Charge and Spin Practionalization Beyond the Luttinger Liquid Paradigm 
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It is well established that at low energies one-dimensional (ID) fermionic systems are described by 
the Luttinger liquid (LL) theory, that predicts phenomena like spin-charge separation, and charge 
fractionalization into chiral modes. Here we show through the time evolution of an electron injected 
into a ID t-J model, obtained with time-dependent density matrix renormalization group, that a 
further fractionalization of both charge and spin takes place beyond the hydrodynamic limit. Its 
dynamics can be understood at the supersymmetric point (J — 2t) in terms of the excitations of 
the Bethe-Ansatz solution. Furthermore we show that fractionalization with similar characteristics 
extends to the whole region corresponding to a repulsive LL. 
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There is a sustained interest in the physics of one- 
dimensional (ID) quantum systems due to recent exper- 
imental advances that allow to access exotic phenomena 
like spin-charge separation and charge fractionalization 
[1]. At low energies these systems are well described 
by the Luttinger Liquid (LL) theory [2] that predicts 
two independent excitations carrying either only charge 
(holons) or only spin (spinous) and propagating with dif- 
ferent velocities, and hence, spin-charge separation. Ex- 
perimental evidences of its existence have been observed 
in quasi- ID organic conductors [3], semiconductor quan- 
tum wires [4] , and quantum chains on semiconductor sur- 
faces [5]. The LL theory also predicts the fractionaliza- 
tion of injected charge into two chiral modes (left- and 
right-going) [6] , a phenomenon recently confirmed exper- 
imentally [7] . Along the experimental advances also the- 
oretical progress was recently achieved pertaining exten- 
sions beyond the LL limit by incorporating nonlinearity 
of the dispersion, leading to qualitative changes in the 
spectral function [8-12] and relaxation processes of ID 
electronic systems [13]. 

Here we show that fractionalization of charge and spin 
beyond the forms described by LL theory takes place 
when a spin- 1/2 fermion is injected into a strongly cor- 
related ID system, namely the t-J model. By study- 
ing the time evolution of the injected wavepacket at dif- 
ferent wavevectors fc, using time-dependent density ma- 
trix renormalization group (t-DMRG) [14-19] different 
regimes were obtained. When k is close to the Fermi 
wavevector kp, the known features from LL theory like 
spin-charge separation and fractionalization of charge 
into two chiral modes result. On increasing k, a fur- 
ther fractionalization of charge and spin appears, in forms 
that depend on the strength of the exchange interaction 
J or the density n. Their dynamics can be understood 
at the supersymmetric (SUSY) point J = 2t ui terms 
of charge and spin excitations of the Bethe-Ansatz solu- 



tion [20-22]. For the region of the phase diagram [23, 24], 
where the ground state corresponds to a repulsive LL, two 
qualitatively different regimes are identified: one regime 
with Vs > Vc and another where < Vc- Here Vc[s) is the 
velocity of the excitations mainly carrying charge (spin). 
For Vs > Vc and k > kp the spin excitation starts to 
carry a fraction of charge that increases with k while Vc 
corresponds to a wavepacket carrying only charge. For 
Vs < Vc and k > kp the situation is reversed and the 
fastest charge excitation carries a fraction of spin that 
increases with k while the wavepacket with Vs carries al- 
most no charge, i.e. in this case spin fractionalizes. 
The Hamiltonian of the ID t-J model is as follows, 

-{-j'^yS^-Si+i-^Uim+ij, (1) 

where the operator cj^ (ci.o-) creates (annihilates) a 
fermion with spin a =t, i on the site i. They are 
not canonical fermionic operators since they act on a re- 
stricted Hilbert space without double occupancy. Si = 
^1 a'^aP^i p is spin operator and rii — ^c^ ^ is the 
density operator. 

We study the time evolution of a wavepacket, corre- 
sponding to a fermion with spin up injected into the 
ground state, by means of t-DMRG [14-19]. The state of 
a gaussian wavepacket \ip) centered at xq, with width A^; 
and average momentum ko, is created by the operator ip^^^ 
applied onto the ground state \G): 

lV')^^jlG)=^^.4lG), (2) 

i 

with 

ip^ = ^e-(a;.-^o)V2A,e*'=«^\ (3) 
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A is fixed by normalization. The time evolved state 
\'<P{t)) by the Hamiltonian (1) determines the spin (s) 
and charge (c) density relative to the ground state as a 
function of time r measured in units of 1/t {H — 1), 

P„(x,,t) = (V^(T)|n,„|^(T)> - (G|n,„|G), (4) 

where a = s,c, riic = n,-^ + Ui^, and riis = - n^. 
Most of the numerical results were carried out on systems 
with L = 160 lattice sites, using 600 DMRG vectors (this 
translates into errors of the order of 10^* in the spin 
and charge density up to times of 50/i) and = 5 
lattice sites (which corresponds to a width ~ 0.067r 
in momentum space). 

We discuss first the time evolution of a wavepacket at 
the SUSY point J = 2t, since here we will be able to iden- 
tify the different portions in which the wavepacket splits 
on the basis of the Bethe-Ansatz solution. Figure 1 shows 



is observed (wavepackets Pi with i = 2,3,4), beyond the 
prediction of the LL theory. 
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FIG. 1. (color online). Time evolution of pc{xi,T) for a 
wavepacket initially at a; = 0, with momentum k = O.Tvr, 
at density n — 0.6, and J — 2t. Charge fractionalizes into 
four wavepackets, one to the left(Pi) and the rest(P2, Pa, Pi) 
to the right. Pi and P3 have the same charge and speed but 
opposite velocities. 

the time evolution of Pc{xi,T) for a density of n = 0.6 
(n = iV/L, with N the number of fermions). The mo- 
mentum of the injected fermion is fc = O.Ttt, i.e. midway 
between kp = O.Stt and the zone boundary. The charge 
(i.e. pc) splits into four fractions, one portion traveling 
to the left and the rest doing so to the right. A splitting 
into chiral modes is expected in the frame of LL theory 
[6], where for an injected right-going fermion, a splitting 

Qq^-* = (1 ± Ka)/2 (where Ka is the so-called LL pa- 
rameter and ('-') corresponds to the right (left) prop- 
agating part) is predicted. The amount of charge (i.e. 
the integral of the wavepacket over its extension) corre- 
sponding to the portion denoted Pi is Qi"'' 0.1 which 
agrees well with the prediction of LL theory, since for the 
parameters in this case, Kc ^ 0.8 [24]. However, at long 
enough times, a further splitting of the right-going charge 
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FIG. 2. (color online). Charge {pc(xi,T), full line) and spin 
{ps{xi,T), dashed Une) densities for J = 2t, n = 0.6, at time 
T = 40, for different values of the momentum of the injected 
fermion. 

Figure 2 displays both Pc{xi, t) (full line) and Ps{xi, t) 
(dashed line) for different values of the initial momentum 
of the injected wavepacket. The arrows indicate the di- 
rection of motion of each packet. As opposed to Pc, Ps 
does not split. In an SU(2) invariant system Kg = 1 
[2], and hence qI"^ = 0, i.e. no left propagating part 
is expected for the spin density. (However, a small de- 
pletion in Ps appears traveling to the left, which would 
correspond to Kg ^ 1. Similar findings were presented re- 
cently [25] and attributed to finite-size effects that require 
exponentially large systems in order to recover Kg = 1.) 
Moreover, part of the charge (P4) is accompanying the 
spin, such that spin-charge separation does not appear 
to be complete. The amount of charge accompanying the 
spin increases as the momentum of the injected fermion 
approaches the zone boundary. These results make al- 
ready evident that injecting a fermion at a finite distance 
from the Fermi energy leads to fractionalization of charge 
beyond the expectations from the LL theory. 

In order to understand the new forms of fractional- 
ization that go beyond the LL frame, we consider the 
excitations corresponding to one-particle addition pro- 
cesses, whose energies can be obtained from the Bethe- 
Ansatz solution [12, 26]. When adding an electron with 
momentum fc, the single particle excitation energy is 
given by ^(fc) = -edqc) - es(9s), where edqc) and 
^sils) are the dispersion relations of the excitations for 
charge and spin, respectively, and the momenta are re- 
lated to the momentum of the incoming particle as fol- 
lows: fc = ±2kp — Qc — qs, where qc € [—qpc,qFc\, and 
qs e [-qps,qps], with qp^ = tt - 2kp and qps = tt - kp 
the pseudo-Fermi momenta for the excitations for charge 
and spin, respectively [22, 27]. 

Figure 3 displays the velocities obtained from t-DMRG 
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FIG. 3. (color online). Full lines: derivatives Va{k) = 
dea(k) / dk of the dispersions obtained by the Bethe-Ansatz 
solution. The symbols correspond to the velocities of the 
different wavepackets identified in Fig. 2: triangles(Pi), 
diamonds(P2), squares(P3) and circles(P4). The orange hor- 
izontal lines stand for the Fermi velocity vpc = dec{<lFc)/dq 
given by the Bethe-Ansatz solution. 

for the different wavepackets (symbols) compared to 
those obtained from Bethe-Ansatz (full lines), as a func- 
tion of the momentum of the injected fermion. The 
velocity of each P.^ is extracted by measuring the po- 
sition of the maximum of the packet at the most con- 
venient time, i.e, at that time where we can resolve Pi 
and the spreading of one packet does not destroy the 
other packets. The wavepackets Pi (triangles) and P3 
(squares) have opposite directions, but the same speed 
and charges ' ~ QcJ': where the charges for the {+) 
branch are labelled by an index corresponding to the re- 
spective wavepackets. The velocity of the wavepacket P4 
(circles) agrees almost perfectly with the one correspond- 
ing to spin excitations. Its determination is best since 
it is the fastest wavepacket, such that it can be easily 
discerned from the rest. The velocity of the remaining 
wavepacket, P2 (diamonds), is more difficult to assess, 
since it overlaps at the beginning with other ones. Nev- 
ertheless, its velocity closely follows the one of charge 
excitations. The wavepackets just described deliver a 
direct visualization of the excitations appearing in the 
Bethe-Ansatz solution, where only two different kind of 
particles are involved: the c and s pseudoparticles with 
their associated bands. The excitation associated with 
spin involves one hole in the c band with fixed momen- 
tum qpc and one hole in the s band with momentum q^, 
where Qs = zt2kp — qpc — k [22]. In fact, the velocities 
of Pi and P3 correspond to the group velocity at both 
pseudo-Fermi momenta iqpc- Furthermore, as shown in 
Fig. 3, the velocity of those fractions is independent of 
the momentum of the injected fermion, in agreement with 
the picture given by Bethe-Ansatz. The dispersion of the 
hole in the s-band gives rise to the velocity displayed by 



the red line in Fig. 3. Similarly, the c line (black line in 
Fig. 3) involves one hole in the s band with fixed momen- 
tum qps and one hole in the c band with momentum qc 
determined in terms of k by Qc — ±2kp — qpg — k. Using 
the same argument as for the s line we can associate the 
P2 packet (diamonds) with the c pseudoparticle. How- 
ever, in this case we cannot observe wavepackets associ- 
ated with spin and velocities corresponding to the group 
velocity at the pseudo-Fermi momenta iqps- We under- 
stand this as due to the fact that = 1, by analogy 
to what we observe in the Kc = 1 case. On the SUSY 
point this case is reached in the limit of vanishing den- 
sity, where the system can be described by a Fermi gas. 
Hence, fractionalization is absent in this limit. 
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FIG. 4. (color online). As in Fig. 3 but for J = 1.75t. The 
fuUines correspond to the SUSY Bethe ansatz. 

Next we depart from the SUSY point and examine how 
fractionalization takes place in the region of the phase di- 
agram where the ground state corresponds to a LL with 
Kc < 1. Figure 4 shows the velocity of the different 
fractions at J = 1.75t, where a slight decrease (increase) 
in the velocity of the spin (density) fraction can be ob- 
served. As shown in Fig. 5, essentially the same features 
are observed as at the SUSY point both for J > 2t and 
J < 2t. In all the cases shown in Fig. 5, where the ve- 
locity of spin excitations («<,) remains higher than that 
of charge excitations (uc) in most parts of the Brillouin 
zone, spin does not fractionalize, as opposed to charge, 
so that the interpretation derived from Bethe-Ansatz re- 
mains valid over an extended region of the phase diagram: 
charge splits into four portions of which one travels with 
the spin wavepacket, and two have the same speed but 
opposite group velocity which does not depend on the 
momentum of the injected fermion. It is tempting to as- 
sign those excitations to states at a pseudo-Fermi surface 
for charge excitations. 

For smaller values of J/t than those in Fig. 5, Vg be- 
comes smaller than Vc- Figure 6 shows that for Vs < Vc 
the role of spin and charge wavepackets experience a 



4 



0.04 
0.03 
0.02 
0.01 


0.04 
0.03 
><, 0.02 
of 0.01 


0.02 
0.01 




J=1.5t,k=0.557l, -c=35 




100 



FIG. 5. (color online). Fractionalized wavepackets for dif- 
ferent values of J/t away from the SUSY point, at a density 
n = 0.6. As in the SUSY case, charge fractionaUzes into four 
pieces, while spin does not, and carries an appreciable amount 
of charge. 



change with respect to fractionahzation. In this case it 
is the spin density that spHts into two fractions, one at- 
tached to the fastest fraction of charge, an another one 
left behind. Again, this new feature is not predicted 
by LL theory. However, no left propagating fraction of 
spin could be observed (excepting the depression due to 
finite-size effects). We therefore expect that this is due to 
SU(2) symmetry and the fact that in this case Kg = 1. 
As shown in Fig. 6, the amount of spin accompanying 
the charge increases as the momentum of the injected 
fermion increases. While such a phenomenon may sug- 
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FIG. 6. (color online). Fractionalized wavepackets a,t J/t — 1, 
n = 0.5, r = 50 and L = 200. In this case, where Vs < fc, 
fractionahzation of the spin density is observed. 

gest as in the lowest panel of Fig. 5 a total recombination 
of charge and spin as the energy increases, it is not total, 
since still a fraction of charge goes to the left, without 
accompanying spin. 

In summary, we have shown through the time evolu- 



tion of an injected spinfull fermion onto the t-J model, 
that charge and spin fractionahzation occurs beyond the 
predictions of the Luttinger liquid theory. A compari- 
son with results from Bethe-Ansatz allowed to identify 
charge and spin excitations that split into components 
at high and low energies. The components at high en- 
ergy reveal the dispersion Cc and es of charge and spin 
excitations, respectively. The components at low energy 
have a velocity that does not depend on the momentum 
of the injected fermion and are very well described by 
states at the pseudo-Fermi momenta of the charge exci- 
tation. This picture can be extended to a wide region in 
the phase diagram of the t-J model as long as the ground 
state corresponds to Kc < 1, and Vs > Vc- In this region 
fractionahzation is observed only in the charge channel. 
However, for Vc > Vs, a region that develops for J/t be- 
low ~ 1.5, the spin density shows fractionahzation. All 
over, the fastest excitation is accompanied by the com- 
plementary one, such that spin-charge separation is for 
them only partial. The other fractions present an almost 
complete spin-charge separation. Finally, we would like 
to remark, that the time evolution helped to visualize 
all the fractions in a more direct way than the study of 
the one-particle spectral function, where previous stud- 
ies [28] only can identify the dispersion corresponding to 
fractions P2 and P4, but not those propagating at the 
pseudo-Fermi points. 
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